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Abstract 

In a natural or artificial channel, the presence of vegetation at the bottom 
and on the wall is usual. The effects produced by the presence of this is 
an alteration of the flow field. The magnitude of the changes in the flow 
field depends on the extension, height, density, and flexibility of the 
vegetation. To analyze the turbulent structures in detail in this work, an 
experimental study was carried out in a laboratory, in a rectangular 
channel with a section of submerged, flexible, and artificial vegetation. 
The measurement region was located upstream of the vegetation zone 
until the end of it, in addition, the asymmetry that the presence of the 
channel wall can induce was determined. The measurements were 
performed with an Acoustic Doppler Velocimeter (ADV) instrument. The 
results allow to identify the development of the mixing layer from the 
average velocity patterns and the Reynolds stress. In the case of the 
average velocity profiles over the vegetation layer, it was possible to 
identify an inflection point where Kelvin-Helmholtz instabilities type were 
identified. In order to analyze the processes in the mixing layer, a 
phenomenological model is proposed based on the analysis of the 
curvature of the averaged velocities. In addition, the contribution to the 
moment transport rate by the identified turbulent structures were 
determined with a quadrant analysis of fluctuating velocities. 
Keywords: Flexible vegetation, pnenomenological model, ADV, mixing 


layer, turbulent structures. 
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Resumen 

En un cauce natural o artificial es usual la presencia de vegetación en el 
fondo y en la pared. Los efectos que produce la presencia de ésta es una 
alteración en el campo de flujo. La magnitud de las alteraciones en el 
campo del flujo depende de la extensión, altura, densidad y flexibilidad 
de la vegetación. Para analizar a detalle las estructuras turbulentas, en 
este trabajo se realizó un estudio experimental en laboratorio en un canal 
de sección rectangular con un tramo de vegetación sumergida, flexible y 
artificial. La región de medición se ubicó aguas arriba de la zona de 
vegetación hasta el final de esta misma, además se determinó la asimetría 
que puede inducir la presencia de la pared del canal. Las mediciones se 
realizaron con equipo de velocimetría acústica de efecto Doppler (ADV). 
Los resultados obtenidos identifican el desarrollo de la capa de mezcla a 
partir de los patrones de velocidad media y los esfuerzos de Reynolds. En 
el caso de los perfiles de velocidad media sobre la capa de vegetación fue 
posible identificar un punto de inflexión donde se identificaron 
inestabilidades tipo Kelvin-Helmholtz. Con el fin de analizar los procesos 
en la capa de mezcla se propone un modelo fenomenológico basado en el 
análisis de la curvatura de las velocidades promediadas. Además, se 
determinó la aportación en el transporte de momento de las estructuras 
turbulentas con un análisis de cuadrantes con las velocidades fluctuantes. 
Palabras clave: vegetación flexible, modelo fenomenológico, ADV, capa 


de mezcla, estructuras turbulentas. 
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Introduction 


In an artificial channel the interaction of the flow with the vegetation 
produces a wake zone and generates regions of turbulence production 
that alter the patterns of average flow, with a remarkably high three- 
dimensional rate. 

In the case of a natural channel the presence of submerged 
vegetation is normal, but changes in its morphology are related to its 
presence and interfere depending on its height, density and flexibility. 
Vegetation generates changes in turbulent structures vertically and 
transversely, as well as important secondary flows. To evaluate the 
magnitude of these changes it is advisable to estimate patterns of average 
flows and Reynolds stresses, as a low frequency characterization 
condition, but also to evaluate the fluctuating velocities statistically. 

This part of the document describes the conditions of submergence, 


bending and density of vegetation. The condition of submergence, 
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according to Palmer (1945) establishes it as the relationship between the 
height of the vegetation with respect to the water depth of the channel. 
The internal composition of this is very important, since it influences in a 
very specific way in the modification of the bottom of the channel. In 
addition, an important parameter to consider, should be the relationship 
between the bending response of vegetation elements in relation to the 
kinetics of the flow field. 

The flexion of vegetation (J) is the dynamic response of the 
structure of the stems of vegetation to the action of the main velocity of 
flow (Okamoto 8 Nezu, 2010; Ghani et al., 2013; Jeon, 2015). This is a 
function of the main velocity and can be classified according to the 
experimental results obtained by Carollo et a/., (2005), and Okamoto 8 


Nezu (2009), this classification is listed below: 


Rigid. The vegetation is maintained in its vertical position and does not 
present oscillations in the top. 

Smooth oscillation. The vegetation has a smooth deformation and has an 
unorganized smooth oscillation. 

Monami. The vegetation has a significant deformation, and a coherent and 
harmonic oscillatory motion is observed. 

Prone. The vegetation has a despondency of all its elements and the 


oscillations are small with high frequency. 


In addition to the bending behavior, vegetation can be classified 
with respect to its location and distribution which can be on the wall or 


bottom of the channel, exclusively or in both places. It's also important 
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to establish the density of vegetation, which can induce different 
behaviors for the same type of relationship between velocity and bending 
(Poggi et al., 2004; Nezu 8 Sanjou, 2008). 

The relationship between vegetation height and water depth, 


defined by Nepf 8 Vivoni (2000) are classified into three types: 


Low flows, the flow passing through the stems and foliage are 
deformed. 

Intermediate flows, have a complete submergence of vegetation, but 
it affects the behavior of the flow field. 

High flows, the flow field is not affected by the presence of vegetation, 
for this condition the vegetation is considered as a uniform bottom 


roughness. 


The study of intermediate flows has become an important topic to 
analyze due the effects that are generated in the mean and turbulent flow, 
since they modify the profiles of average velocity generating coherent 
turbulent structures near the top of the vegetation and these vortices are 
mainly large-scale, generating an exchange zone of mass and 
momentum, between the vegetation top and outer zone, that is, on the 
vegetation, a condition similar to a boundary layer. In engineering 
practice submerged vegetation is analyzed as a condition of uniform 
resistance, similar to a high flow. Then, to evaluate its effect, the usual 
concept of the resistance of the flow by effect of the submerged 


vegetation and its effects on the roughness parameters is applied, 
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obtaining general values of Manning's roughness coefficient as a function 
of the submergence ratio (Járvelá, 2004; Wilson, 2007). 

In this work we will carry out a more localized study of the effects 
of submerged vegetation, considering that the development of the flow 
field on the top of the vegetation is similar to the condition of two flow 
fields at different velocities and with a layer of interaction between these 
and that produces small perturbations, which induce an instability of the 
flow of the Kelvin-Helmoholtz (KH) type, so above large-scale vortices are 
generated, as Finnigan (2000) demonstrate. To analyze in greater detail 
the generation of turbulent structures in principle it can be said that have 
a boundary layer and this develops longitudinally, but Ghisalberti 8 Nepf 
(2002) demonstrate that flows on submerged vegetation can be analyzed 
as a mixing layer, rather than as a boundary layer. The velocity profile of 
a mixing layer approximates the shape of a hyperbolic tangent profile (Ho 


8: Huerre, 1984), as shown in Figure 1. 
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Figure 1. Velocity type profile in a blend layer (Ghisalberti € Nepf, 
2002). 


According to Poggi et al. (2004) in the vegetation zone three regions 
can be identified, according to the size of the vortices, as shown in Figure 
2. The first region is mainly dominated by small-scale vortices (von 
Kármán vortices). The second region, which is near the top of vegetation 
is dominated by Kelvin-Helmholtz vortices generated by instability at the 
inflection point. The third region is similar to a boundary layer-like 


development zone. 
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Figure 2. Three vortices scales in a channel with submerged vegetation 
(Poggi et al., 2004). 


Sukhodolov € Sukhdolova (2006) performed field measurements to 
verify the applicability of the analogy of the mixing layer in a natural 
channel, on a patch of Sagittaria sagittifolia implanted with a certain 
arrangement varying their density. The results confirm the effectiveness 
of the mixing layer analogy. Sukhodolov € Sukhdolova (2006), and 
Maltese et al. (2007) focused on the spatial pattern of the turbulent 
structure that develops over submerged vegetation. In addition, they 
applied quadrant analysis, revealing that the dominant structure is the 
expulsion on the vegetation top. 

Ghisalberti € Nepf (2002) revealed that flows with submerged 


vegetation can be modeled ¿as a mixing layer, rather than a boundary 
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layer zone, because the velocity profiles that occur in a channel with 
submerged vegetation approximate a hyperbolic tangent profile, see 
Figure 1. In addition, from the average velocity profiles they identify an 
inflection point, which implies that the flow ¡is susceptible to Kelvin- 
Helmholtz instabilities. Nikora et al. (2008) examine the effect of 
submerged vegetation on hydraulic resistance in a variety of vegetation 
patch patterns. Okamoto et al (2012) performed continuous inkjet 
experiments to evaluate the structure of mass transport in an open 
channel flow with rigid vegetation models, and a rigid strip configuration 
varying density. Obtaining the effects of coherent vortices on vertical 
turbulent diffusion. Okamoto € Nezu (2013) examined the transition from 
a boundary layer upstream the vegetation patch to a mixing layer type 
flow developed over the vegetation and developed a phenomenological 
model for the development process in the flow with submerged 
vegetation. In these experiments they used the PIV (Particle Image 
Velocimetry) technique, and the vegetation elements were composed of 
rigid strips. 

As mentioned above, the characteristics of mean flow and 
turbulence in open channel flows with submerged vegetation have 
received a lot of attention. However, detailed information on the evolution 
of the mixing layer and the coherent turbulent structures generated by 
the presence of vegetation is not yet available. For this reason, the 
present study considers the effect of a patch of flexible vegetation of finite 


length and an intermediate flow to identify turbulent structures, as well 
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as their longitudinal and transverse evolution from a phenomenological 
model. The analysis zone starts upstream of the patch to the patch 
boundary, but cross-sectional measurements were made and vertically 


the area below the vegetation top to the outer flow was considered. 


Materials and methods 


The experimental station is a channel with fixed slope S, of rectangular 
section, with a length of 10 m, a width of 0.60 m and with a height of 0.60 m. 
The walls of the channel are made of glass and steel bottom. The 
representation of the vegetation elements is of a flexible plastic material. 


The hydraulic conditions of the experimental study are shown in Table 1. 
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Table 1. Hydraulic conditions. Where Q is the flow rate, U = Q/A is the 
average velocity of the flow, Sy is the slope of the channel, Re = UR,,/v is 
the Reynolds number, Fr = U/,/gR,, is the Froude number, R, = 4/P is 
the hydraulic radius, A is the area of the cross-section, P is the 
perimeter of the cross-section, v is the kinematic viscosity and g is the 


acceleration of gravity. 


h(cm) | h,(cm) | h/h, |Q(1s*) U(cms”) So Re Fr 


33 11 3 41.20 20 0.005 32698 0.11 


The reference system for sampling from the experimental station is 
shown in Figure 3. Where the coordinate system is the triad of values of 
xi, for ¿=1,2,3, and is named as longitudinal, vertical and transverse 
coordinate respectively. The domain of the measuring space in meters is 
0(x¿) € [-0.5, 3.7] x [0.0, 0.33] x [0.0, 0.3]. 
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Figure 3. Experimental station. 


For sampling of the instantaneous velocities in the three 
components (u;), an Acoustic Doppler velocimetry (ADV) device was used 
at a sampling frequency rate of 100 Hz during a period of T=60s. The 
ADV device is of the Nortek6) or also called Vectrino Profiler”" (Vectrino 
11). For the processing of the experimental data, we considered the 
Reynolds decomposition model u; =U,+u;, considering the temporal 
averaged velocities in each component as U, and the fluctuating velocities 
as u;. 

The vegetation configuration is a patch with a length of 3.70 m and 
the arrangement of flexible elements is linear, with a separation between 


elements in longitudinal direction is l,,=0.025m, and in transverse 
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direction is l,, = 0.025 m, as shown in Figure 4 and Figure 5, respectively. 


The start of the vegetation patch is considered in x, = 0. 


X2 
Flow 


Smooth bed zone Vegetation zone 


h =0.33 m 


—=—-).50 SN o X1 


Figure 4. Vegetation zone. 


X3 


X1 


Figure 5. Patron of vegetation elements. 
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The sampling points were from the center of the channel to the right 
bank in five sections, indicated by the 4,B,C,D,E letters, see Figure 6 and 
longitudinally 43 points were analyzed and 19 vertically, see Figure 7a. 


The measurement points in the vertical are shown in Figure 7b. 


Lo oosr o. ....... . » 
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Figure 6. Measurements in the transverse direction. 
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Figure 7. Measurement points with the ADV equipment, with a total of 


4085 points. a) Sampling points in the plane (x;,,x3) with a separation of 
0.10 m in the direction x, and 0.06 m in the direction x3; b) Sampling 
points in the plane (x,,x>2) with a separation of 0.015 m in the direction 
x2. In the vegetation zone there are 7 points inside and 12 above the 


vegetation. 


According to Vargas-Luna et al. (2015), the actual vegetation can 
be simulated with cylindrical bars. With this, it was proposed to use 
cylindrical bars to simulate vegetation with the following physical 
characteristics: a height of h, = 0.11 m and a diameter of 0 = 7.3 mm, the 
material of the vegetation elements is silicone. 

For the characterization of the vegetation were used the parameters 
of vegetation density, Af, flexural rigidity, J, and the ratio of 
submergence, h/h,. The dimensionless vegetation density, Af, was 


determined with the expression proposed by Okamoto € Nezu (2013): 
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_ y Af ¡hp 
y == (1) 


Where 4; is the total frontal area of the vegetation elements, 4), is 
the area occupied by the vegetation. In the case of flexural rigidity, J, Of 


vegetation elements, the equation proposed by Wilson et al (2003), was 
used, and ¡s: 


J=El=2== (2) 


Where E is the modulus of elasticity, ! the second moment of the 
area (also called the moment of inertia). On the right side of the 
expression, we have that F is the relationship between the force, 6, ¡it is 
the deflection of the element, L it is the length of the vegetation element 
(Figure 8). Then, we get a value of J = 3.42x10"* Nm* and a value 2; = 


1.37. 
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Figure 8. Diagram of the vegetation element with deflection. 


For this study we simulated a type of real vegetation called 
Lespedeza sericea with a length of h, = 10— 20cm and with a J = 5x107”* 
1x1073 Nm? (Kouwen 8 Li, 1980; Nepf 8 Vivoni, 2000). The parameters 
used in this experiment, such as, the submergence ratio, the density, the 
flexural rigidity and the pattern of the vegetation elements, was based on 
various laboratory experiments carried out by various authors (Palmer, 
1945; Wilson et a/., 2003; Nezu € Sanjou, 2008; Okamoto 81 Nezu, 2013). 
Linear and stepped patterns are the most used in the laboratory, because 
they have a better representation of the effect of vegetation in a channel 


with intermediate flow (Vargas-Luna et al., 2015). 


Tecnología y ciencias del agua, ISSN 2007-2422, 
14(2), 27-82. DOI: 10.24850/j-tyca-14-02-02 


del Agua. Open Access bajo la licencia CC BY-NC-SA 4.0 
(https://creativecommons.org/licenses/by-nc-sa/4.0/) 


a LN) Check for updates 
OPEN ACCESS 
Tecnología y 


CienciaszAgua 


Results 


According to the above, the flow condition being analyzed is an 
intermediate flow. In addition, for a better understanding of the results 
we considered a dimensionless space of the sample domain 0(¿€, = x,/h,) € 
[-4.54, 33.63] x [0.0, 3.0] x [0.0, 2.72], the top of the vegetation is held for ¿, = 
1, and the start of the vegetation for ¿, = 0. For the cross-sectional 
sections shown in Table 2. 


Table 2. Values of ¿z for the measurement points in the transverse 


direction. 


Section A B C D E 


Eg =x3/hy 272 2.18 1.63 1.09 0.54 
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Mean flow characteristics 


Of the measurements obtained, the fields of mean flows were first 
analyzed. Figure 9 shows the temporal mean velocity contours in the main 
direction U,. The mean velocity at the sampling points (¿;) was determined 


with the following expression. 


UE =5f (E mdn (3) 
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Eg =2.72 Eg =2.18 


U,(m/s) 


Figure 9. Contours of mean velocity in the main direction U,. Where the 


dashed line represents the top of the vegetation ¿, = 1. 


Due to the presence of the vegetation, a deceleration of the main 
velocity is observed near the vegetation top ¿, = 1, which is compensated 
by an increase in this in the outer zone. This distortion of the flow field 
causes an inflection point right in the vegetation top ¿, = 1, the above 


generates a Kelvin-Helmholtz instability and therefore vortices of different 
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scales are presented, this observation resembles results presented by 
Nepf 8 Vivoni (2000). 

In addition, Figure 10 shows the mean velocity profiles for all cross 
sections, observing a hyperbolic profile. These flow field characteristics 
are similar to those presented in terrestrial vegetation, presented by 
Finnigan (2000). The instability is maintained in the transverse direction, 


but the thickness of the mixing layer decreases due to the presence of 
the wall. 


E =-454  £,=0.90 E, = 6.36 E =1181  £ =17.27 £=22.72  f,=28.18 E, = 33.63 


Cross — Section 

—e— E, = 1,63 
—e—¿, = 2.72 —+— f, = 1,09 
—e— f, = 2,18 —+— E, = 0,54 


U, mio 


Figure 10. Mean velocity profiles in the main direction, U¿. Where the 


dashed line indicates the top of the vegetation. 
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The mean velocity profile in ¿, = —4.54, see Figure 10, in the zone 
near the free surface (¿€, > 2.5), an acceleration of the flow is observed by 


the backwater effect caused by the vegetation. 


Reynolds stress 


Analyzing Reynolds stress —uju, of the experiments (Figure 11), the 
development of the mixing layer can be observed, which is very 
representative for the central area (¿3 = 2.72), but this is decreasing as the 
influence of the walls is had, é€z < 2.72. 
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Figure 11. Reynolds stress contours —u;¡u; . Where the dashed line 


represents the top of the vegetation ¿, = 1. 


Nepf € Vivoni (2000), Nezu € Sanjou (2008) define the size of the 
mixing layer from the Reynolds stress profiles —u; us, where the upper and 


lower limit are defined by the 10 % of the maximum stress |-ujuz]___, that 


is, as shown in the following expression: 
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Where y, is the upper limit and y, is the lower limit of the mixing 
layer. The upper and lower limits of the mixing layer for ¿¿ = 2.72 and é, = 
0.54 are shown in Reynolds stress profiles, according to the profile of 
Reynolds stresses for the longitudinal position, and in the cross section 
¿1 =33.63you have a mixture layer thickness of and ¿3 = 2.72for ty, = 
0.19 mof ¿ =0.54. tm, = 0.15 m (Figure 12). With this it is observed that at 
the center of the channel you have a larger mixture layer size than you 


have glued to the wall. 
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Figure 12. Reynolds stress profiles —u¡u). Where the blue and magenta 
dashed lines represent the longitudinal spatial variation of the thickness 
of the mixing layer in the center section ¿3 = 2.72 and near to the wall 
¿3 = 0.54, respectively, and the black dashed line represents the top of 


the vegetation ¿, = 1. 
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Mixing layer 


According to Raupach et al., Finnigan and Brunet (1996), and Ghisalberti 
8: Nepf (2002), the velocity profile of a mixing layer can be analyzed with 


the following expression: 


2" tanh (2) (5) 


For better understanding, Figure 1 shows the variables in the mixing 
layer. Where U,,, = U,(y,) + U,(y,)/2, AU, is the velocity difference for each 
cross section: AU, = U¿(y,) — U¿(y,), U¡(y,) and U;, (y, ) indicate the velocity 
averaged at the upper and lower limit of the mixing layer respectively, y, 
and y, is the distance from the bottom of the channel to the value of the 
Reynolds stress with zero value below and above the maximum Reynolds 
stress value respectively, y = (y, + y,)/2, this analysis is developed for 
each velocity profile in the longitudinal and cross-sectional sections. The 


momentum thickness 6 is defined with the following expression: 


oo za) Jas (S 
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Figure 13 makes a comparison between the theoretical model to 
evaluate the hyperbolic tangent velocity profiles of a pure mixing layer 


and the experimental data. 
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Figure 13. Comparison of the mean velocity profiles with the hyperbolic 
tangent profile of a mixing layer for each cross section for different 


positions in the longitudinal direction. 


According to the results it is observed that the velocity profiles 
obtained in the experimentation resemble a mixing layer profile, 


therefore, it is concluded that a mixing layer is formed in an intermediate 


flow. 


Phenomenological model 


To evaluate the evolution of the mixing layer longitudinally for each 
analysis section, it is proposed to apply a division of the flow field into 


four zones are described in Table 3 and shown in Figure 14. 
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Table 3. Division of the phenomenological model zones, where X, is the 


beginning of analysis (zone upstream of the vegetation), Xp, is the 


beginning of the divergent zone, Xp is the beginning of the 


development zone, Xzp is the beginning of the fully developed zone, and 


X¡ is the end of the field of analysis. 


Zone Flow condition Analysis space 
I Initial [X;, Xp ] 
II Divergent [XpL, Xp] 
III Development [Xp Xzp] 
IV Fully developed zp, XA] 
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Figure 14. Phenomenological model of flow development through a 
zone of submerged vegetation (Okamoto € Nezu, 2013; Huai € Qian, 
2015). 


Upstream of the vegetation zone (¿, < Xp,) the velocity profiles 
resemble a profile found in an open channel flow, since the vegetation 
does not affect the flow. At the beginning of the vegetation (¿, = 0) there 
is a high resistance to the flow and initiates an increase in the velocity 
above the vegetation and a deceleration within it. That is, part of the flow 


is diverted at the beginning of the vegetation zone. Zong 8 Nepf (2010) 
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call this phenomenon as divergent flow. The end of the divergent flow is 


denoted as Xp. 


Downstream the divergent flow zone (¿, > X»p), the velocity profiles 
have a high velocity gradient, as shown in Figure 1. This gradient creates 
an inflection point at €, =1, inducing a Kelvin-Helmholtz instability, 
generating downstream vortices. The growth of these vortices is seen in 
the Reynolds stress profiles. This zone is defined as a development zone, 
the end of this is defined as Xzp. 

Downstream of the development zone (¿, > Xp) the size of the 
vortices remains constant. This area ¡is called a fully developed zone. 

Figure 15 shows the effect of vegetation on the main velocity 
averaged vertically U,, in the zone without and with vegetation. To 


determine this averaged velocity, the following expression is used: 


Ub == 45 "UE E, E = ctedag (7) 
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Figure 15. Variation of the mean velocity in the main direction, 


upstream and over vegetation. 


Observing in the zone upstream of vegetation, —5< é¿, < 0, there is 
a velocity with little variability, 9U,,, /dx = 0. For the central zone of the 
channel ¿; = 1.63, 2.18, 2.72, the flow velocity increases to a maximum value 
in a range 0<gé,<5, downstream this maximum value the velocity 
decreases, because the mixing layer reaches the free surface, as shown 


in Figure 11, being a zone of high turbulence. However, for the zone near 
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to the wall ¿¿=0.54,1.09, the increase in velocity remains constant, 
because the mixing layer does not reach the free surface, therefore, the 
turbulence zone is smaller and the deceleration of the flow is not 
generated, see Figure 11. 

The starting point of the divergent flow, X;,, and the length of the 
Divergent Zone (Zone 11), Lrp = Xpp — Xp,, is determined from the vertical 
velocity averaged in the vertical, U,, . Applying equation (7) evaluates the 
divergent flow, as shown in Figure 16. Observing a rapid increase in 
vertical velocity in -2.73 <¿, <0, where the maximum value is at the 
beginning of the vegetation zone (¿, = 0) for each section. The start of the 
divergent flow, X,, is determined when the vertical velocity increases, 
that is, when dU,, /dx, > 0, therefore, for the five sections X,, = —2.73. 
However, the start of the development zone, Xp, for each section varies 
as it approaches the wall. X,, is obtained up to the point where the vertical 
velocity approaches zero, U,, =0, as shown in Figure 16. Then Table 4 


shows the length of Zone IT, Lyp. 


Tecnología y ciencias del agua, ISSN 2007-2422, 
14(2), 27-82. DOI: 10.24850/j-tyca-14-02-02 


del Agua. Open Access bajo la licencia CC BY-NC-SA 4.0 
(https://creativecommons.org/licenses/by-nc-sa/4.0/) 


0.05 


0.04 


0.03 


-0.01 


-0.02 


Tecnología y 


OPEN Que: 


CienciaszAgua 


1) Check for updates 


XrDE 


XFDp XrDc XFDp XFDA 


—0— ¿3 = 2.72 
—4— ¿q = 2.18 


A — Eq = 1.09 
—P— Ez = 0.54 


Sas 


V 


SH ¿3=1.63 | 


5 


10 


15 


ó1 


25 


30 


Figure 16. Variation of the average velocity vertically for each section. 


Table 4. Divergent flow length, Lp. 


Zonevés 2.72 2.18 1.63 1.09 0.54 
XrL 2.73 2.73 2.73 2.73 2.73 
Xrp 13.63 12.72 10 7.27 2.72 
Lrp 16.36 15.45 12.73 10 5.45 


To better visualize the Diverging Flow Zone, Figure 17 shows the 
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zone with the highest vertical velocity U,. 
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Figure 17. Contours of mean velocity in the vertical direction U,. Where 


the dashed line represents the top of the vegetation ¿, = 1. 


Downstream of the divergent flow zone the velocities vertically 
become small, observing negative velocity values, U, < 0. Then for ¿, > 
Xp the mixing layer it grows gradually until the size of the mixing layer, 


tm, remains constant, that is, when dt,,/dx, = 0, is located in the start of 
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the fully developed zone, Xzp. Figure 12 shows the spatial variation of the 
mixing layer. Then for each section, Xzp is obtained with the spatial 


variation of the size of the mixing layer, as shown in Figure 18. 
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Figure 18. Spatial variation of the mixing layer thickness, tn, for each 


section. 
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For comparison, Figure 19 shows the variation of momentum 
thickness, 06, which is a measure of growth of the shear layer. Noting that 


tr and 6 grow rapidly downstream from the start of the vegetation zone. 
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Figure 19. Spatial variation of the momentum thickness, 6, for each 


section. 
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Then the length of the Development Zone is obtained with Lp = 


Xzp — Xfp. Table 5 shows the Zone III lengths for each section. 


Table 5. Length of development zone, L;p. 


Zonevéz 2.72 2.18 1.63 1.09 0.54 
Xp 13.63 12.72 10 7.27 2.72 
Xzp 19.09 20.91 21.82 25.45 27.27 
Lzp 5.45 8.19 11.82 18.18 24,55 


Finally, downstream of Zone III, €, > Xzp, is the Fully Developed 
Zone (Zone IV), L¿p, to where the vegetation zone ends. Figure 20 shows 
the phenomenological model of flow development for this test, which 


shows the variations of the four zones in the transverse direction, xz. 
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Figure 20. Phenomenological model. Where the black dashed line 


indicates the start of the vegetation zone, ¿, = 0. 


It is observed that the length of Zone III, L,», for the center of the 
chamnel, ¿3 = 2.72, is less than the zone near to the wall, ¿¿ = 0.54. This 
indicates that in the center of the channel there is a greater resistance, 
since in the center there is a deceleration of the flow is high compared to 


the zone near to the wall, this is observed in Figure 15. 
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Quadrant analysis 


A quadrant analysis was performed for instantaneous Reynolds stress 
4 (0)us(t) for the analysis of turbulent structures that contribute to 
momentum transport. The basic principle of quadrant analysis is to 
decompose Reynolds stresses into four types of flows, called Q,, Q,, (and 
Q, into the plane u; —u; as a function of the signs of fluctuation of the 


main u; and vertical us velocity as shown in Figure 21. 


Figure 21. Quadrant diagram u; — u». 
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Reynolds stress quadrant RS, is defined as follows (Okamoto et al., 
2012; Okamoto 8 Nezu, 2013): 


z Trrs 
RS; = lim 7 f, (uz)! de (8) 


Where T is the duration of the measurement. If (u;,uz) exist in 
quadrant Q;, then /; = 1, otherwise I; = 0. 
Each quadrant (u;,u,) corresponds to the following coherent events: 
Q, (u; > 0,us > 0) : External interaction, which is the upward movement of 
the high-speed flow. 

Q, (u¿ <0,us > 0) : Ejection, which is the upward movement of the low- 
speed flow. 

Q3 (u¿ > 0,u5 < 0) : Internal interaction, which is the downward movement 
of the low-speed flow. 

Q, (u; <0,u5 < 0) : Sweep, which is the downward movement of the high- 
speed flow. 

Figure 22 shows the vertical distribution of the Reynolds stress 
quadrants, RS;, for the four zones of the phenomenological model in ¿¿ = 
2.72. In addition, the plane of the fluctuating velocities (u; — us) for ¿, = 
1.2, where they indicate the percentage of contribution of each quadrant. 
In the initial flow zone, é¿, = -2.72, it is observed that the values are small 
and the percentage of contribution of turbulent structures to the 
momentum transport is similar, that is, no organized structures are 


observed. 
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Figure 22. Reynolds stress quadrants for the central section of the 
channel ¿3 = 2.72. a) initial zone (¿, = -2.72); b) divergent zone 
(€, = 5.45); Cc) developed area (¿, = 17.27); d) fully developed zone 


($1 = 30). 
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In the Divergent Flow Zone, é, = 5.45, over the vegetation zone 
(€, > 1), it is observed that the dominant structure is the sweep, RS,, see 
Figure 22.b. 

In the Development zone, é, = 17.27, the dominant structures are 
the expulsion and sweep by effect of the vortices in the mixing layer, 
which are the low-speed upward and high-speed downward movements 
respectively Huai et al. (2019). This is observed in £, = 1.5, where the 
maximum value of Reynolds stresses is found (Figure 22c). 

In the Fully Developed Flow zone, é, = 30, the ejection and sweep 
structure increase approximately 11 % with respect to the Development 
zone. These results indicate that the contribution of momentum transport 
between the vegetation zone and the zone above it is mainly due to the 
conditions of expulsion and sweeping. 

To visualize the turbulent structures throughout the flow region, 
Figure 23 shows the contours of the Reynolds stress quadrants for all 
sections. It is observed that for all sections the structures that contribute 
most to the transport of momentum is the structure of ejection and 
sweep. In addition, it is observed that the maximum values are in the 


mixing layer zone. 
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Figure 23. Contours of the Reynolds stress quadrants, RS;, for each 


cross section ¿z, evaluated over space xy — Xz. 


Once the turbulent structures with the greatest contribution have 
been found, Figure 24 compares the contributions of the sweep and 


ejection movement, using the Equation (9): 
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Figure 24. Relationship between sweep and ejection structures. 
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(9) 


Applying Equation (9) for each section we obtain the relations 
between sweep and ejection Rs,¿, observing a relationship Rs,g > 1 for all 


cross sections, therefore, the structure with the greatest contribution to 
vertical momentum transport is due to sweep motion. 
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Conclusions 


In this study, turbulence measurements were performed in a channel with 
a patch of vegetation of finite length to examine the creation and 
evolution downstream of turbulent structures in the mixing layer above 
vegetation. 

The mean velocity profiles confirm the presence of an inflection 
point caused by Kelvin-Helmholtz type instability similar to theoretical 
analysis when two potential flows of different velocity are presented. By 
vertically distributing Reynolds stresses, the mixing layer was identified, 
as well as the variation in size along the vegetation patch and transverse 
direction. Noting that the size of the mixing layer stops growing in the 
Fully Developed Zone (x, > Xzp), because the mixing layer extends to the 
free surface. In addition, the size of the mixing layer for ¿z = 2.72 (section 
A) is approximately 20 % larger compared to ¿z = 0.54 (section E), this is 
clearly seen in the contours of Reynolds stresses, see Figure 11 because 
the high-speed flow field is larger than the low-speed, as shown in Figure 
9. 

With the phenomenological model of the flow, the four zones of the 


flow were identified for the five analysis sections, obtaining the lengths of 
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these zones for the sections, as shown in Figure 20. Looking for ¿z¿ = 0.54 
(section E) requires a greater distance for the mixing layer to fully develop 
compared to ¿z = 2.72 (section A). This is visualized in the magnitude of 
Reynolds stresses. 

A quadrant analysis was applied to identify the turbulent structures 
that contribute to the momentum transport, it is visualized that for the 
Initial Flow Zone an organized structure is not observed. In the 
Development and Fully Developed Zone, ejection and sweep structures 
appear near the top of the vegetation. These results indicate that the 
contribution of momentum transport between the vegetation zone and 
the zone above it is mainly due to the ejection and sweep structures. 

With the methodology used in this paper it is possible to identify the 
presence of a mixing layer and the characteristics of the turbulent 
structures present in it. In addition, with the use of a phenomenological 
model applied to several sections, a map was generated that shows the 
variation of the four zones along a channel with vegetation, which can be 
used for channels with vegetation with characteristics similar to those 
used in this experiment, such as the submergence ratio, flexural rigidity, 


dimensionless density and the main mean velocity. 
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